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Abstract 



We develop the Inverse Scattering Transform (1ST) method for the Degasperis- 
Procesi equation. The spectral problem is an sl(3) Zakharov-Shabat problem with 
constant boundary conditions and finite reduction group. The basic aspects of the 
1ST such as the construction of fundamental analytic solutions, the formulation of 
a Riemann-Hilbert problem, and the implementation of the dressing method are 
presented. 



1 Introduction 

The Degasperis-Procesi (DP) equation 



03 ■ u t - utxx + 3ku x + iuu x - 3u x u xx - uuxxx = 0, (DP) 

where k > is a constant, was first discovered in [11] in a search for asymptotically inte- 
grable PDEs. Equation (DP) is a bi-Hamiltonian system and admits interesting traveling 
wave solutions [12]. It arises as a model equation in the study of two-dimensional water 
waves propagating over a flat bed [39, 10, 36]. Given the intricate structure of the full 
governing equations for water waves, it is natural to seek, in various physical regimes, 
simpler approximate model equations. In the case of two-dimensional waves propagat- 
ing mainly in one direction, two parameters appear in the non-dimensionalisation of the 
governing equations: the wave-amplitude parameter e (indicating how close the waves 
are to a fiat-surface flow) and the long-wave parameter 8 (measuring the ratio of the 
approximate wavelength to the average water depth) [38]. The relative sizes of these 
two fundamental parameters determine the different physical regimes for water waves. 
The most studied regime is the shallow-water regime (also called the small amplitude or 
long- wave regime) for which S <C 1 and e ~ S 2 . In this parameter range, due to a balance 
between nonlinearity and dispersion, various integrable systems like the Korteweg-de 
Vries (KdV) equation arise as approximations to the governing equations [2]. The inte- 
grability of these equations implies that the powerful method of inverse scattering can be 
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used to obtain detailed qualitative (and even quantitative) information about the wave 
dynamics. In particular, since linear water wave theory cannot provide an approximation 
of solitary waves (see the discussion in [8] ) , an important outcome of studies of the KdV 
equation was a deeper understanding of the dynamics of solitary water waves [38] . 

The shallow amplitude regime is, however, not appropriate for the study of large 
amplitude waves, whose behavior is more nonlinear than dispersive. To model such 
waves, which are characterized by a relatively large value of e, it is natural to investigate 
the parameter regime in which J < 1 and e <~ 5. Also in this parameter range is a 
reduction to a simple wave equation at leading order possible, but since the dimcnsionless 
parameter e is larger, the nonlinear effects are stronger than in the shallow water regime. 
Of particular interest in this regard is that a stronger nonlinearity could allow for the 
occurence of wave-breaking — a fundamental phenomenon in the theory of water waves 
that is not captured by the KdV equation (see the discussion in [7]). In the shallow- 
water, moderate amplitude regime [8 <C 1, £ ~ S] several model equations can be derived 
as approximations to the governing equations for water waves. However, among the 
various equations that arise in this way, there are only two which admit a bi-Hamiltonian 
structure (see [10, 35]): the Camassa-Holm (CH) equation [4] and the DP equation (DP). 
For these equations, u(x,t) represents the horizontal velocity of the water at a certain 
depth at time t and position x; the depth ratios for CH and DP are different [39, 10] 
and are encoded in (DP) in the positive non-dimensional parameter n > 0. In the limit 
k — ^ the solitary waves of (DP) are peakons and are given explicitly by 

u c (x,t) = ce~ lx ~ ctl , x,teR, 

where c > denotes the speed of the wave [12, 44]. Due to the presence of a peak at 
the wave crest, these waves have to be understood as weak solutions [42]. This feature 
is characteristic for waves of great height (i.e. waves of largest amplitude that are exact 
solutions of the governing equations for water waves, see the discussion in [6, 8, 51]). The 
dynamics of the peakon interactions for DP was elucidated in [45, 46, 13]. 

In this paper, we will develop an inverse scattering approach for smooth localized 
solutions to (DP). More precisely, we consider solutions u(x, t) of class C 1 in time and 
of Schwartz class regularity with respect to the spartial x-variable (i.e. the solution is 
smooth and decays to zero faster than any polynomial as |x| — > oo). Moreover, we will 
assume that the solution satisfies the following inequality initially (at time t = 0): 

q = u — u xx + k > 0. (1.1) 

Note that well-posedness and global existence of solutions for (DP) holds within the class 
of Schwartz functions if the initial data satisfy (1.1), and in this case the validity of (1.1) 
is ensured at any later time t > (see [14, 31, 43]). 

The paper is organized as follows. In Section 2 we present the Lax pair formulation of 
(DP) and study the symmetry properties of the isospectral problem and of the scatter- 
ing matrix. Section 3 is devoted to the associated Riemann-Hilbert problem, while the 
Zakharov-Shabat dressing method is implemented in Section 4. 



2 Spectral problem 
2.1 Lax pair 

Equation (DP) admits the Lax pair formulation [12, 33] 

- q( 3 <P = 0, 




(2.2) 

+ utp x - u x (p = 0, 

where £ £ C is the spectral parameter and ip(x, t) is a scalar function. Third-order 
spectral problems appear as eigenvalue problems for the Boussinesq, the Sawada-Kotera 
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and the Kaup-Kuperschmidt equation [41] and they were first investigated by Kaup 
[40] who established the analyticity properties of the fundamental solutions, constructed 
scattering data and presented an inverse scattering approach in the cases when only 
the continuous spectrum is present and when only the bound state spectrum is present. 
The third order problem (2.2) however is a weighted spectral problem which requires 
additional care. For example, a weighted second order spectral problem arises in the 
inverse scattering for the CH equation, its analysis is described in details in [9] . 

Our aim is to develop the inverse scattering approach for (DP). The Lax represen- 
tation (2.2) however is inconvenient for this purpose. It turns out that the problem 
simplifies a lot if the Lax representation is written in the form of Zakharov-Shabat (ZS) 
type spectral problem [54, 50, 49, 48, 15, 25, 17, 18, 19]. This allows us to take advantage 
of the existing inbuilt symmetries of the considered equation. For this reason we write 
the above Lax pair in matrix form as 

f** =l *> (2.3) 

where 

/_! ( Q\ ( u + u x + ^s -C«-£ \ 

L= C , M=\ -W £"C« 

\(q 1/ y -quC'+q^-Ux+u^ q+Ux+Uxx _„ _ ^ + ^ 

and (p(x,t) is a SL(3) - matrix-valued function, whose columns, considered as vectors, 
represent the three linearly independent solutions of the matrix equation. Let G(x, t) be 
a SX(3) matrix and let us take a change of variables <p = Gip. It transforms (2.3) into 

(ip x = Lip, 
\ip t = Mip, 

where 

L = G- 1 LG-G~ 1 Gt, M = G- 1 MG-G~ 1 Gt. 



Letting uj = e 2 " 1 / 3 and 



-1/3 
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1 





1 




UJ 


u? 












UJ 


= Q* 


(1 


3ql 







'uj 0\ 1 / uj + 1 1 

uj 2 and Q* = -(l-w) 1 oj + 1 \ . (2.4) 

v 1/ 3 \uj + 1 1 

Let us change the variables according to 



y = x + 



1 



da;' 



= fq{x)\ 
dx V k J 



q(x)\ 1/3 

K 



(2.5) 



The t-variable can be viewed as an additional parameter rather than a second indepen- 
dent variable. This will be clear in the following section and is due to the fact that 
the ^-dependence of the scattering data can be explicitly computed in relatively simple 
form. For the sake of simplicity in what follows we usually omit the t-dependence of 
the variables, unless this dependence is necessary for the computations. The spectral 
problem 

ipx + {Q- (q 1/3 J)ip = (2.6) 
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can be written in the form 



4> y + {Q- AJ)^ - 0, (2.7) 

where 

and lim 2; _j.± 00 h(x) = 1. 

Suppose that x = X{y). It is possible to recover g(x) from h(X(y)). First, we notice 
that asymptotically y — > x when x — > — oo. Since 



f [(- 



-g(x)V/ 3 _ 1 



dx 



is an integral of motion [12], for x — > oo we have that x and y differ only by a constant. 
Let us introduce for convenience 



Then /, and therefore q, can be recovered from h(y) = h(X(y)) by solving the first order 
differential equation that follows from (2.8): 

+ %)/ = !, /(±oo) = l. (2.10) 
We obtain the solution in parametric form 

q(X(y)) = K f(y), x = X(y) = y + J* (-^ - l)d W ', (2.11) 
or, with the convolution interpreted in the sense of distributions [34], 

/oo 
S(x-X(y))f 2 (y)dy, (2.12) 
-oo 

where the function X(y) is defined in (2.11) and f(y) is the solution of (2.10). Finally, 
h(y), which is actually h(X(y)), can be obtained from the scattering data of the spectral 
problem (2.7) which is a ZS-type spectral problem, however with constant boundary 
conditions, since Q(y) — Q*h(y) and lim y _,.±oo h(y) = 1. For other ZS spectral problems 
for multicomponent systems with constant boundary conditions we refer to [3, 23, 15]. 



2.2 Automorphisms 

The specific form of the potential Q in (2.7) is due to the symmetry of the problem 
under the action of three distinct automorphisms. In other words, the fact that Q is 
determined by a single real (scalar) function, rather than 6 complex functions (which 
is the case for an arbitrary sl(3) potential) is a consequence of its invariance under 
one Z3 automorphism and two Z2 automorphisms. The automorphisms lead to the 
reduction of the independent components of the potential and their action extends to 
the spectrum and the eigenfunctions. They form a group, known as a reduction group 
[47, 24, 25, 22, 20, 21, 27, 28, 29, 30]. 
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2.2.1 Z3 automorphism 

The spectral problem (2.7) has a manifest Z 3 symmetry: 



CL^XjC- 1 = L{\), (2.13) 

where 

C = 




Indeed, one can check that 

CJC- 1 = — J, CQC- 1 = Q 
u 

from where (2.13) follows immediately. Furthermore, one can verify that 

CM(u\)C- 1 = M(X), (2.14) 

so that the Z 3 automorphism (2.13) is an automorphism of the graded Kac-Moody algebra 
where L(A) and M(A) take their values [16]. This holds for all automorphisms of the 
spectral problem. 

2.2.2 Z 2 automorphisms 

The spectral problem possess two additional Z 2 automorphisms, one of which reflects the 
reality of h(y). The first one is 



TL(wA)r- 1 = L(A), (2.15) 

where 




The second Z 2 automorphism is 

ZL^-'^Z- 1 = -L(X), (2.16) 

where 





>\ 



























and the dagger stands for a matrix Hermitian conjugation. 

One can easily check that all Z 2 and Z3 automorphisms are compatible, in the sense 
that the order of their application to the potential does not matter. 

3 Construction of Fundamental Analytic Solutions 

3.1 Asymptotic behavior 

Let us define the asymptotic values 

L oc (X)= lim L(y,X), M QO (A) = lim M(y, A). 

y— !-±oo y^±oo 
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Then 

Loo = At/ ~ Q 
where Q* is defined in (2.4). We also find that 



3iU 2 X 2 LU 2 {CO - 1) A - LO LU 2 {1-LU)\-L 

M oc (A) = ^ ( (1 - u)X-lu 2 3wA 2 (w-l)A-w ) (3.17) 

tj(u-l)X-u cj (1 - oj) A - w 2 3A 2 



3A 3 

and that 



[Loo, Moo] =0. 



Since and M M commute, they can be simultaneously diagonalized. Let U(X) be a 
SL(3) matrix such that 

ioo(A) = U(X)A(\)U- 1 (\), MUX) = U(X)A(X)U- 1 {X), 

with 

A(A) = diag(A 1 (A),A 2 (A),A 3 (A)), A(\) = diag(A 1 (\),A 2 (\),A 3 (\)), 

where Ai, A 2 , A 3 and A\,A 2 , As are the eigenvalues of L M and Mqo respectively. 

The eigenvalues of Loo (A) = A J — Q* arc the solutions A(A) of the characteristic 
equation 

A 3 - A - A 3 = 0. 



(3.18) 



Introducing a new spectral parameter k such that 
we have 

A 3 = 3- 3 / 2 (fc 3 + -l) (3.19) 
and the following expression for the eigenvalues of L^: 

Aj(k) = T 1/2 (io 1 k + -L-y (3.20) 

Furthermore, X(k) has the property X(ojk) = ojX(k) and also 

X(k) -> 3~ 1/2 k when |fc| -> oo. (3.21) 

The characteristic polynomial of the matrix 3k _1 A 3 Moo, cf. (3.17), is 

P( w ) = w 3 -3w- 27A 6 + 2, (3.22) 

with roots 



Wj(k) = u j k 2 + u j k 2 , 
where A(fc) is given in (3.18). Thus the eigenvalues of M m are 

A,m - 

It remains to determine the ordering of the eigenvalues Aj that is consistent with the 
ordering of Aj , the eigenvalues of . To this end we consider the asymptotic expressions 
when k — > oo. Then 



kJ , r \/3k t9 



G 



or 



A,^, Aj^^uP. (3.24) 



Thus, from (3.23) we have 



A, W . ^^±^T1; (3,5) 

3.2 Scattering matrix 

Consider the modified Lax pair 

U V = I4, L = \J-Q(y) 

The compatibility condition holds for any choice of matrix replacing M 00 (X), i.e. the 
modified Lax pair gives rise to the same equation as the original Lax pair. For convenience 
we use the spectral parameter k. Let ^(y, t, k) be the solutions of (3.26) such that 

lim ^(y^k) = U(k)e A W y U-\k). 

y— »±oo 

Note that due to the modified second equation in (3.26) the asymptotic values do not 
depend on t. The two solutions ip + {y,t, k) and tp~(y,t,k) are not linearly independent, 
i.e. they are related by a linear transformation when k is a point of the spectrum. Thus, 
the expression tp + (y 7 t, k)tj)~ (y, t, k) depends on t and k, but not on y. (From now on we 
write B for the inverse of a matrix B). We define the scattering matrix T(t, k) by 

T(t,k) = U(k)i>+(y,t,k)^-(y,t,k)U(k), (3.27) 

The explicit form of the i-dependence of T is quite simple: 

Proposition 3.1 The time- evolution of the scattering matrix is given by 

T(t, k) = e A(k)t T(0, k)e- A{k)t . (3.28) 

Proof. From (3.27) it follows that 

^+UT = ^~U. (3.29) 

Differentiating both sides with respect to t, we obtain 

ip+UT + i{j + UT t = iptU. 

Replacing ip t from the t-part of the Lax pair (3.26), we find 

(M^+ - V + Moo)C/T + ijj+UTt = (Mip~ - ifr-Mo^U. 

In view of (3.29) this becomes 

-V+MooC/T + 4> + UT t = -t^+UTUM^U. 

We conclude that T(t, k) evolves according to T t — — [T, UMooU], or 

T t = [A,T]. 

Therefore the time-evolution of the scattering matrix is given by (3.28). □ 
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Figure 1: Domains of analyticity and integration contours. 



3.3 Fundamental analytic solutions 

An important role in the theory of inverse scattering play the so called fundamental 
analytic solutions (FAS) of the spectral problem. We will explain the construction of 
FAS for the system (2.7) or rather for the related spectral problem 

Sv + QS + [S,Hk)]=0, Q = U(Q-Q*)U, (3.30) 

in the domains v = 1, . . . , 6 occupying the space in the complex fc-plane separated 
by the rays {l v : arg(fc) = (y — l)f } with v = 1, 6 (see Figure 1). It is easy to notice 
that the components of Q(y) are functions in the Schwartz class. 

Let us introduce an ordering for each sector fi„ as follows. We will say that 



r<s iff 5ft[A r (fc)] < 5ft[A s (fc)] when k g Q v . (3.31) 

One can verify that 1 < 2 < 3 since for k g fii, SR[Ai(fc)] < SR[A 2 (fc)] < SR[A 3 (fc)]. Similarly, 
1 < 3 < 2, etc. Then we can prove the following result: 

2 2 

Proposition 3.2 The solution £ v (y, k) of the following system of integral equations 



CMV = 5 rs - r e-^- A ^y-y'^Q(y\k)C(y',k)] rs dy', r>s; (3.32) 

J-oo v 

C s (y,k) = - [\-^- A ^y-y\Q{y',k)^{y'Mrsdy l , r<s. (3.33) 

Joe v 

is also a solution of spectral problem (3.30) and is analytic for k e £L V . 

Proof. The first part of the statement follows by direct computation. The analyticity 
follows from the fact that the real part of A r (fc) — A s (k) is nonnegative when k g VL V and 
r>s etc. □ 

V 

By direct computation we prove the following useful result. 
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Proposition 3.3 The solution £(y,k) of (3.30) is related to a solution ip(y,k) of (2.7) 
via 

£(y t k) = U{k)i;(y,k)U(k)B(k)e- A ^ y 
where B(k) is some nondegenerate matrix. 

Using the two propositions we compute the asymptotics 

lim [e-W> v C(y, k)e A{k)v ] rs = 0, r < s, fc G ft„. 

y— 7-oc v 

Therefore 

lim e- A ^y^(y, k)e A(k ^ = T~{k)D+(k), k G Q„, (3.34) 

where T~(k) is a lower -triangular matrix with units on the diagonal with respect to the 
^-ordering and D+(k) is a diagonal matrix. Similarly, 

lim e -A(fc)s/^( y) fc) e A ( fe )f =S+(k), ken u , (3.35) 

y — > — oo 

where S+(k) is an upper-triangular matrix with units on the diagonal with respect to 
the ^-ordering. 

A more-subtle result, the proof of which follows exactly the lines of the proof of 
Theorem 3.3. from [26], gives the asymptotics for k G l v : 

Proposition 3.4 When k G l v , both [S+{k)] rs = and [T-(k)] ra = if $t(A r (k)) ^ 
»(A a (fc)). 

In other words, on the line k G l v nonzero could be only the entries that have the 
numbers r , s of the two eigenvalues that define l v as 5i(A ro (A;)) = 5i(A So (fc)). For 
example, when v = 2, r = 2, ,s = 3. Notice that if in r < s 0} in the neighboring 

v — 1 

fi„ (l v separates fi„_i and Q v ) the ordering changes such that r > s . Therefore, if the 

V 

matrix S^_ 1 (k) is upper triangular with respect to the v— \ ordering, it is lower-triangular 
with respect to the ^-ordering and we will denote it as S~(k), etc. 

The asymptotics on the ray k G l v therefore can be written in the form: 



lim e- A(k ^C(y,k)e A(k ^ = T~{k)D+{k), (3.36) 

y—>oo 

lim e- A(k ^C(y,k)e A(k ^ = S+(k), (3.37) 

y—y — oo 

lim e- A ^r _1 (y,fc)e A(fe)2/ = T+{k)D-{k), (3.38) 

J/^OO 

lim e- A ^r _1 (y,fc)e A(fe)!/ - S~(k), (3.39) 



where the matrices S^, (resp. S~ , T~ ) arc upper-triangular (resp. lower triangular) 
with units on the diagonal (with respect to the ^-ordering) and the matrices D ± are 
diagonal. They provide the Gauss decomposition of the scattering matrix with respect 
to the ^-ordering [16, 26], i.e. 

T v (k)=T-(k)Dt(k)S+(k)=T+(k)D-(k)S-(k), kel v (3.40) 

where T~,S~ are lower-triangular matrices with units along the diagonal, T + ,S + are 
upper-triangular matrices with units along the diagonal, and are diagonal matrices. 
From the definition (3.27) of T and (3.40) it follows that on the lines k G l v , 

i>+(y,k)U(k)T-(k)D+(k) = ^-(y,k)U(k)S+(k), (3.41) 
i) + {y,k)U{k)T+{k)D-{k) = i>-(y,k)U(k)S-(k). (3.42) 
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Now, using the results from Proposition 3.3 and (3.36)- (3.42) we can relate the eigen- 
functions with tp^ for k e l v : 



C(y,k) = U{k)i,-{y,k)U{k)S+{k)e-^y 

= tj(k)i; + (y,k)U(k)T-(k)D+(k)e- A ^ v (3.43) 

?-\y,k) = U(k)4,-(y,k)U(k)S-(k)e- A ^y 

= U(k)i; + (y,k)U(k)T+(k)D-(k)e- Mk)y . (3.44) 

The symmetries (2.13), (2.15) and (2.16) impose constraints to the eigenfunctions, 
the scattering matrix and its factors e.g. 

CC +2 {y,ojk)C = £ v (y,k), (3.45) 

rf-"(i,,wfc)f - r(y,fc), (3.46) 

Z[?- v (y,-k)]*Z = C(y,k), (3.47) 

etc. where the ^-indices should be taken modulo 6. 

Thus, independent is only the data on one of the lines, say lo, and all the rest can be 
recovered from (3.45) - (3.47). 

4 Riemann-Hilbert problem 

Any two solutions £ v (y, k; t) and £"~ 1 (y, k; t) are related on k e l v due to (3.43) - (3.44): 

r(y,M) = C-\y,k;t)G v {y,k;t), (4.1) 
G v {y,k;t) = e Mk)y+Aik)t S-(k)S+(k)e~ Mk)y - A{k)t k e l v , (4.2) 
lim C(y,k;t) = i. (4.3) 

k— >oo 

We can show that the relations (4.1) - (4.3) constitute a Riemann-Hilbert problem 
(RHP) for the matrix-valued functions £"{y,k;t), each being analytic for k € Q u . In 
other words, the solution of the original Inverse Scattering Problem reduces to a RHP 

Proposition 4.1 The RHP (4-1) - (4-3) has a unique solution ^ u {y,k;t), with an ana- 
lytic continuation for k € 

For simplicity we consider the case where £ v (y,k;t) do not have any singularities in 
the fc-planc. Suppressing y and i-dependence for convenience where possible, we can 
write an analytic continuation for ^ v (y, k) (k e Sl v ) as follows: 

i/=i 



on,, 



where the orientations of the contours <9f2„ are shown on Figure 1 and C v belong to the 
infinite circle. Due to (4.3) we have 

' C(y,k')dk' f C(y,k')dk' , , +1 f ldk' x 



2ni ^ W (i/ k' — k Ji„ +1 k' — k } Cv k' — k j 
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where the orientation of l v is always from to oo. From (4.1) it follows that 

If k approaches l n from the left and right domains correspondingly, with the Sohotski- 
Plcmelj-type formulae we obtain that 

PW . , + it^wp-W -i] + JL g w . /_ Mi) , 
- . - ^- wi^w -n + ^ (t «v. X r ~' ( T- ( r """"' j ■ 

fc g J„ t/= 1,2. ..,6 (4.7) 
From (4.7) and (4.1) we finally obtain a system of integral equations for ^ v ~ 1 (y, k): 

- . + ^(e^./ ^ 1( ~'-^ ). 

fc e Z„ »y = l,2...,6. (4.8) 

The solutions of (4.8) provide ^ r '~ 1 (y,k) when k £ l v and then the analytic continu- 
ation is given by (4.6). Apparently the scattering data are provided by G u (y, fc; t) given 
in (4.2), i.e. by S~(k) and <S+(fc), fc € l v . Using the Z 2 and Z 3 automorphisms (2.13) 
- (2.16) i.e. (3.45) - (3.47) we can restrict the scattering data to 5f(fc) and 5^(fc), 
fc S Zi. However, according to Proposition 3.4, Si~(fe) has only one nontrivial component 
[S^ (fc)]i2, which is exactly the only nontrivial component of S± (k). Again, the automor- 
phisms (3.45) - (3.47) relate S^ik) to i5j~(fc) which restricts the minimal set of scattering 
data to only one function, say [S' 1 h (fc)]i2. This function is the analogue of the reflection 
coefficient used for example in the Inverse Scattering for the KdV equation [48] . 

The potential Q(y,t) of the scattering problem can be recovered from the following 
result: 

Proposition 4.2 If £ v (y,h;t) are the solutions of the RHP (4-1) - (4-3) then 

X "(y,k;t) = U(k)C(y,k;t)e A ^y (4.9) 

satisfy (2.7) with 

Q{y, t) = - lim X(k) ( X v (y, k; t)J X v (y, k; t) - J) . (4.10) 
Proof. With arguments similar to those given in [49, 50, 16] one can prove that 

Q(y, t) = - lim X(k) U"(y, k; t)j£"(y, k; t) - j) . (4.11) 
From (4.11) and the definition of Q in (3.30) it follows that 

Q(y, t)=Q*- lim X(k) ( X v (y, k; t)J?(y, k; t) - U(k)JU(k)) . (4.12) 
Next, we notice that 

lim (A(jfe) - X(k)J) = 0, 

k^-oc 
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giving 

lim (u(k)A(k)U(k) - X(k)U(k)JU(k)) = 0, 
lim (l x - X(k)U(k)JU(kj) = 0, 

or 

lim (\(k)J - Q* - X(k)U(k)JU(k)) = 0, 

i.e. 

Q* = lim X(k) (j - U(k)JU(kj) . (4.13) 

Clearly (4.10) follows from (4.12) and (4.13). □ 

Corollary From (4-10) and the fact that tr{{Q*) 2 ] = 2, one can find 

Hy,t)= l -tr{QQ*) (4.14) 

and then f(y,t) can be computed from (2.10) and q(x,t) from (1.1). Finally, u(x,t) can 
be obtained from (2.11) and (2.12). 

5 Zakharov-Shabat dressing method 

The Zakharov-Shabat dressing method [49, 50, 53, 54] allows the explicit construction 
of a solution with singularities x"(y, fc; t), starting from a given regular solution of the 
RHP, say X S(y,M): 

X v {y, fc ; t) = g(y, k; t)\l{y, k; t). (5.1) 

Note that in our case the analyticity regions for x" do not coincide with those for due 
to the nontrivial [/(fc)-factors in (4.9). 

The dressing factor g is analytic in the entire complex plane, with the exception of the 
points of the discrete spectrum. We make the following assumptions in our construction 
of a dressing factor. First, we allow only simple poles of g and g. For simplicity, instead 
of the fc-dependence we are going to revert to the A - dependence, having in mind the 
previously defined relation fc(A). The automorphisms act on g(X) as on a group element: 



Cg(y,LoX)C = g(y,X), (5.2) 
Tg(y,uX)t = g(y,\), (5.3) 
Zgt(y,-X)Z = g(y,X), (5.4) 

From these symmetries it follows that if g or g have a pole at, say, Ao, then they have 
also poles at — Ao, ±wA , ±w 2 Ao, ±Ao, ±wAo, ±cj 2 Ao. It also follows that det(g) = 1. 
Our next assumption will be that A can be chosen real, so that the following choice of 
g is possible 

ff(y> A;t)=I + £^M, (5.5) 

j= i i 

* Za\(y,t)Z 

g(y,X;t)=l-J2 / . (5-6) 
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for some residues aj(y,t) where 



Aj = w i+1 A , Hj = -Xj. (5.7) 

The property (4.3) of the dressed solution is preserved. Following the ZS construction 
from [54] (see also [16, 5]), we represent the residues in the form 

3 

a P = ^\ n i)Rw( m p\ ( 5 - 8 ) 

R jp = ^ (5.9) 
where \rij) is a vector-column, and (rrij\ is a vector-row. The last is defined as 

( mj \ = (moj\x U Hv,*j)- (5-10) 

where (m j | is a constant vector and x^ is the solution of the adjoint problem, analytic 
at Xj. We define 

|nj> ee Z\rrij) = X 0) Oj)Kj), where \nrij) = ((m i |) T (5.11) 
and X^K^j) — Z[x^\^j)V %i Le - the constant vectors are related as follows: 



\n 0j ) = Z\m 0j ). (5.12) 

X^ is an eigenfunction, analytic at fij, since it is obtained by the automorphism (2.16) 
from x^\ analytic at Xj (and Aj and fj,j are related by the same automorphism). With 
these definitions one can easily check that the matrix R (5.9) is Hermitian: R = B) . 
The residues of g in (5.6) can be computed and represented in the form 

3 

Zal(y,t)Z = J2K)R Pj (^j\ (5-13) 

One can now verify that gg =1. Indeed, this is satisfied iff at any singular point, say 
A = [i v the corresponding residues satisfy 



3 

a, 



3 (Za]z) 

3- ' 



etc. The identities (5.14) can be verified by (5.8) - (5.13). Our construction for g (5.5) - 
(5.6) satisfies also the automorphism (5.4). The automorphism (5.2) gives an additional 
relation between the constant vectors \n j). 



\n ,j+i) = C\n 0j ). (5.15) 
From (5.3) we have the further restriction 



r|n i> = |n 02 ), r|n 03 ) = |n 03 ). (5.16) 
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The relations (5.15) show that only one of the vectors, say |noi), determines the others 
(and also (moj\ due to (5.12)). The components nni-j (j = 1,2,3) of the vector |no,i) 
are not independent: as a consequence of (5.16) and (5.15) they satisfy noi ; i = fioi-2 
and rioi ; 3 = noi;3- Thus, if noi ; 3 ^ we can take noi ; 3 = 1 and then only one complex 
number, n 01;1 = p determines 

|"o,i) = (Po,Po,l) T (5.17) 

and therefore |no,j) an d \mo,j)- 

Let us denote Qo = Q*ho, where h is the 'undressed' potential. Then we have the 
following equation for g: 

g y + Qg- gQo + \[g,J] = 0, (5.18) 

satisfied identically for any z, i.e. any A. This equation is satisfied for the construction 
(5.5) due to the fact that 

(( m j\)v - (mj\(Q - XjJ) = 0, (5.19) 
(\nj))v + (Qo-HjJ)\nj) = 0. (5.20) 

The identity (5.18) for A — > oo gives 

Q = Qo + [J,ai+a 2 + a 3 ]. (5.21) 

Multiplication of (5.21) by Q* followed by taking of a trace gives 

1 3 

h(y,t) = h (y,t) + - £ pj K|[Q*,J]K>. (5.22) 
p,j=i 

The reality of (5.22) can be checked once again by using the introduced properties of 
R, \rij) and 

Thus, to each discrete eigenvalue, which represents an 'action' variable and is deter- 
mined by only one real value, Ao (5.7) one can put into correspondence a conjugated 
'angle' variable that is given by the independent component of the associated constant 
vector, po (5.17). This accounts for the scattering data related to the discrete spectrum. 
Since A is real, in the k - plane the corresponding discrete spectrum value ko is also 
on the real line. It seems that ko is on the continuous spectrum, since l\ and 1 4 are on 
the real line. However, the continuous spectrum for the actual spectral problem (2.7) 
does not consist of the rays l v shown in Figure 1, but rather of the lines l v through the 
origin, where arg(fc) = 7r/6 + ^7r/3. This is due to the relation between the corresponding 
eigenfunctions (4.9) involving nontrivial fc-dependence in the factor e A ^ v . 

To repeat the dressing procedure N times we will need TV copies of the same type 
scattering data, {po.j, Aoj} j — 1, . . . , N. The result will be the iV-soliton solution. 

For example, the one-soliton solution can be obtained as follows. One can start the 
dressing from the 'trivial' solution u = 0, i.e. h = 1. Then there is a global analytic 
solution of the spectral problem, symbolically e L ^( k )y+ M ^( k ) t i i. e . 

*(y,fc;i) = U{k)e K ^y +A{k ^U(k). 
Since the last factor does not depend on y and t we can take simply (for all sectors) 

Xo (y,k;t) = U(k)e A ^y +A ^ t 1 
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The entries of U(k) are 



U lp (k) 
U 2p (k) 
U 3p (k) 



(uj 2 - 1)A + (uj - l)A p + cj 2 , 
(uj - 1)A + (uj 2 - l)A p + w, 
3(A 2 + \A p + A 2 ) - 1. 



The potential for the one-soliton solution is given by (5.22), then f(y,t) (t is viewed 
as an additional parameter rather than a second independent variable) can be computed 
as a solution of linear first order ODE (2.10). Then q(x 7 1) from (1.1) and therefore u(x, t) 
can be obtained from (2.11) and (2.12). 

6 Discussion 

In the presented analysis we formulated the inverse scattering problem for the DP equa- 
tion. We defined a set of scattering data for the problem: 

(i) on the continuous spectrum the coefficient [S^(k)]i2, k G l\; 

(ii) on the discrete spectrum the set {p j, \ j}, j = 1, . . . , N, where N is the number 
of the discrete eigenvalues. 

The scattering data uniquely define the potential h(y,t) and therefore a general N- 
soliton Schwartz-class solution u(x,t) of DP equation. As it is the case for all systems, 
integrable by the Inverse Scattering Method, the mapping between the solution and the 
scattering data allows the interpretation of a generalised Fourier transform [1, 25, 26, 52]. 

The peakon solutions (peaked solitons) appear in the limit n — > [44] , although such 
limit in the space of the scattering data will require further considerations. 

The iV-soliton solution of the Degasperis-Procesi equation is obtain in [45] by Hirota's 
method. E.g. the 1-soliton solution is 



u(y,t) 



x 




(6.23) 



(6.24) 



where the quantities £j are 





(6.25) 



where and y i0 are constants, representing the scattering data, 




as well as 71 = y (2ai+i)(ai-i) are a ^ so constants depending on the scattering data. 
The quantities, related to the 1-soliton solution (6.23) can also be computed: 




h(y,t) = 



cosh 2^i + 



^-cosh£ 
01 s 



£1 A. 1/1 sinh£i 



cosh£i + 



2+ul 



y (cosh&+£) 
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One can establish that z/j that determines the dispersion relation (6.25) is related to 
the spectral parameters, 



v, 



) 



where fcj is a the real discrete eigenvalue, corresponding to the real Ao,i = \{k{). The 
computation of the soliton solution by dressing method requires development of additional 
techniques addressing the technical difficulties arising in the computation. 

There are interesting multidimensional versions of the DP equation. They are in 
general nonintegrable, but admit singular (peakon-type) solutions [32] . 

Acknowledgements 

A.C. and R.I. acknowledge the support of the G. Gustafsson Foundation for Research in 
Natural Sciences and Medicine (Sweden). R.I. is also supported by Science Foundation 
Ireland, Grant 09/RFP/MTH2144. J.L. acknowledges support from a Marie Curie Intra- 
European Fellowship. The authors are grateful to Prof. V. S. Gcrdjikov and Dr. G. Gra- 
hovski for many valuable discussions. 

References 

[1] M. J. Ablowitz, D. J. Kaup, A. C. Newell, and H. Segur, The inverse scattering transform 
— Fourier analysis for nonlinear problems, Stud. Appl. Math. 53 (1974), 249 -315. 

[2] B. Alvarez-Samaniego and D. Lannes, Large time existence for 3D water-waves and asymp- 
totics, Invent. Math. 171 (2008), 485-541. 

[3] V. Atanasov and V. Gerdjikov, On the multi-component nonlinear Schrodinger equation 
with constant boundary conditions, in Gravity, Astrophysics, and Strings 05, eds. P. P. 
Fiziev and M. D. Todorov, St. Kliment Ohridski University Press, Sofia, 2006. 

[4] R. Camassa and D. Holm, An integrable shallow water equation with peaked solitons, Phys. 
Rev. Lett. 71 (1993), 1661-1664. 

[5] J. L. Ciesliriski, Algebraic construction of the Darboux matrix revisited, J. Phys. A: Math. 
Theor. 42 (2009) 404003. 

[6] A. Constantin, The trajectories of particles in Stokes waves, Invent. Math. 166 (2006), 
523-535. 

[7] A. Constantin and J. Escher, Wave breaking for nonlinear nonlocal shallow water equations, 
Acta Mathematica 181 (1998), 229-243. 

[8] A. Constantin and J. Escher, Particle trajectories in solitary water waves, Bull. Amer. 
Math. Soc. 44 (2007), 423-431. 

[9] A. Constantin, V. Gerdjikov and R. Ivanov, Inverse scattering transform for the Camassa- 
Holm equation Inverse Problems 22 (2006), 2197 - 2207. 

[10] A. Constantin and D. Lannes, The hydrodynamical relevance of the Camassa-Holm and 
Degasperis-Procesi equations, Arch. Ration. Mech. Anal. 192 (2009), 165-186. 

[11] A. Degasperis and M. Procesi, Asymptotic integrability, in Symmetry and Perturbation 
Theory, edited by A. Degasperis and G. Gaeta, World Scientific (1999), pp. 23-37. 

[12] A. Degasperis, D. Holm and A. Hone, A new integrable equation with peakon solutions, 
Theor. Math. Phys. 133 (2002), 1461-1472. 

[13] A. Degasperis, D. Holm and A. Hone, Integrable and non-integrable equations with peakons, 
in Nonlinear Physics: Theory and Experiment (eds: M. Boiti et al.) World Scientific Pub- 



lishing 2007, 37 - 43. 



16 



[14] J. Escher, Y. Liu and Z. Yin, Global weak solutions and blow-up structure for the 
Degasperis-Procesi equation, J. Fund. Anal. 241 (2006), 457-485. 

[15] L. D. Faddeev and L. A. Takhtadjan, Hamiltonian methods in the theory of solitons 
(Springer Verlag, Berlin, 1987). 

[16] V. Gerdjikov, Algebraic and analytic aspects of soliton type equations, Contemp. Math. 
301 (2002), 35-68. 

[17] V. Gerdjikov, Basic aspects of soliton theory, in Sixth International Conference on Geom- 
etry, Integrability and Quantization, June 2004, Varna, Bulgaria, ed: I. Mladenov and A. 
Hirshfeld, SOFTEX, Sofia 2005, pp. 78-125. 

[18] V. Gerdjikov, Generalised Fourier transforms for the soliton equations. Gauge covariant 
formulation., Inverse Problems 2 (1986), 51-74. 

[19] V. S. Gerdjikov, The Zakharov-Shabat dressing method and the representation theory of 
the semisimple Lie algebras, Phys. Lett. A 126 (1987), 184-188. 

[20] V. S. Gerdjikov, G. G. Grahovski and N. A. Kostov, Reductions of A-wave interactions 
related to low rank Lie algebras I: Z 2 - reductions, J. Phys. A 34 (2001), 9425-9461. 

[21] V. S. Gerdjikov and G. G. Grahovski, On N-wave type and NLS type systems and their 
gauge equivalent: generating operators and the gauge group action, Proc. NAS Ukraine 50 
(2004), 388-395. 

[22] V. S. Gerdjikov, G. G. Grahovski and N. A. Kostov, On the multi-component NLS type 
equations on symmetric spaces and their reductions, Theor. Math. Phys. 144 (2005), 1147- 
1156. 

[23] V. Gerdjikov and P. P. Kulish, The multicomponent nonlinear Schrodinger equation in 
the case of nonzero boundary conditions, Zap. Nauchn. Sem. Leningrad. Otdel. Mat. Inst. 
Steklov 131 (1983), 34-46 (in Russian); J. Math. Sci. 30 (1985), 2261-2269 (in English). 

[24] V. Gerdjikov, N. Kostov and T. Valchev, Soliton equations with deep reductions. General- 
ized Fourier transforms, in Topics in contemporary differential geometry, complex analysis 
and mathematical physics, pp. 85-96, World Sci. Publ., Hackensack, NJ, 2007. 

[25] V. Gerdjikov, G. Vilasi and A. Yanovski, Integrable Hamiltonian Hierarchies. Spectral and 
Geometric Methods, Lecture Notes in Physics 748, Springer, Berlin - Heidelberg, 2008. 

[26] V. Gerdjikov and A. B. Yanovski, Completeness of the eigenfunctions for the Caudrey- 
Beals-Coifman system, J. Math. Phys. 35 (1994), 3687-3725. 

[27] G. G. Grahovski and M. Condon, On the Caudrey-Beals-Coifman System and the Gauge 
Group Action, J. Nonlm. Math. Phys. 15 (2008), Suppl. 3, 197-208. 

[28] G. G. Grahovski, On the reductions and scattering data for the CBC system, in Geometry, 
Integrability and Quantization III, Eds: I. Mladenov and G. Naber, Coral Press, Sofia, 
2002, pp. 262-277. 

[29] G. G. Grahovski, On the reductions and scattering data for the generalized Zakharov- 
Shabat systems, in Nonlinear Physics: Theory and Experiment. II, Eds: M. J. Ablowitz, 
M. Boiti, F. Pempinelli and B. Prinari, World Scientific, Singapore, 2003, pp. 71-78. 

[30] G. G. Grahovski, V. S. Gerdjikov, N. A. Kostov, V. A. Atanasov, New integrable multi- 
component NLS type equations on symmetric spaces: Z4 and reductions, in Geometry, 
Integrability and Quantization VII, Eds: I. Mladenov and M. De Leon, Softex, Sofia (2006), 
154-175. 

[31] D. Henry, Persistence properties for a family of nonlinear partial differential equations, 
Nonlinear Anal. 70 (2009), 1565-1573. 

[32] D. Holm, M. Staley, Wave structure and nonlinear balances in a family of evolutionary 
PDEs, SI AM J. Appl. Dyn. Syst. 2 (2003), 323-380. 



17 



[33] A. N. W. Hone and J. P. Wang, Prolongation algebras and Hamiltonian operators for 
peakon equations, Inverse Problems 19 (2003), 129-145. 

[34] L. Hormander, The analysis of linear partial differential operators. I. Distribution theory 
and Fourier analysis, Springer- Verlag, Berlin, 2003. 

[35] R. I. Ivanov, On the integrability of a class of nonlinear dispersive wave equations, J. 
Nonlinear Math. Phys. 12 (2005), 462-468. 

[36] R. I. Ivanov, Water waves and integrability, Philos. Trans. Roy. Soc. London A 365 (2007), 
2267-2280. 

[37] R. I. Ivanov, On the dressing method for the generalised Zakharov-Shabat system, Nucl. 
Phys. B 694 (2004), 509-524. 

[38] R. S. Johnson, A modern introduction to the mathematical theory of water waves, Cam- 
bridge University Press, Cambridge, 1997. 

[39] R. S. Johnson, Camassa-Holm, Korteweg-de Vries and related models for water waves, J. 
Fluid Mech. 455 (2002), 63-82. 

[40] D. J. Kaup, On the inverse scattering problem for cubic eigenvalue problems of the class 
tpxxx + GQipx + 6R1P = Xip, Stud. Appl. Math. 62 (1980), 189-216. 

[41] D. J. Kaup, The legacy of the 1ST, in The legacy of the inverse scattering transform in 
applied mathematics (South Hadley, MA, 2001), pp. 1-14, Contemp. Math., 301, Amer. 
Math. Soc, Providence, RI, 2002. 

[42] J. Lenells, Traveling wave solutions of the Degasperis-Procesi equation, J. Math. Anal. 
Appl. 306 (2005), 72-82. 

[43] Y. Liu and Z. Yin, Global existence and blow-up phenomena for the Degasperis-Procesi 
equation, Comm. Math. Phys. 267 (2006), 801-820. 

[44] H. Lundmark and J. Szmigielski, Multi-peakon solutions of the Degasperis-Procesi equation, 
Inverse Problems 19 (2003), 1241-1245. 

[45] Y. Matsuno, The 7V-soliton solution of the Degasperis-Procesi equation, Inverse Problems 
21 (2005), 2085-2101. 

[46] Y. Matsuno, Multisoliton solutions of the Degasperis-Procesi equation and their peakon 
limit, Inverse Problems 21 (2005), 1553-1570. 

[47] A. V. Mikhailov, The reduction problem and the inverse scattering method, Physica 3D 
(1981), 73. 

[48] S. P. Novikov, S. V. Manakov, L. P. Pitaevskii, and V. E. Zakharov, Theory of solitons: 
the inverse scattering method, New York: Plenum, 1984. 

[49] A. B. Shabat, The inverse scattering problem for a system of differential equations, 
Funkcional. Anal, i Prilozen. 9 (1975), 75-78 (in Russian). 

[50] A. B. Shabat, An inverse scattering problem, Diff. Uravneniya 15 (1979), 1824-1834 (in 
Russian) . 

[51] J. F. Toland, Stokes waves, Topol. Methods Nonlinear Anal. 7 (1996), 1-48. 

[52] T. Valchev, On the Kaup-Kupershmidt equation. Completeness relations for the squared 
solutions, in Ninth International Conference on Geometry, Integrability and Quantization, 
June 8-13 2007, Varna, Bulgaria, ed: I. Mladenov and M. de Leon, SOFTEX, Sofia 2008, 
pp. 1-12. 

[53] V. Zakharov and A. Shabat, A scheme for integrating the nonlinear equations of mathe- 
matical physics by the method of the inverse scattering problem - I, Func. Anal. Appl. 8 
(1974), 226-235. 

[54] V. Zakharov and A. Shabat, Integration of nonlinear equations of mathematical physics 
by the method of inverse scattering - II, Func. Anal. Appl. 13 (1979), 166-174 (English 
translation) . 



18 



